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^y^ Ab initio calculations including electron correlation are still extremely costly except for the smallest 
'««.^toms and molecules. Therefore, our purpose in the present study is to employ a bond-order cor- 
' telation approach to obtain, via equilibrium molecular energies, molecular dissociation energies and 
l/~])heats of formation for some 20 molecules containing C, H, and O atoms, with a maximum number 
C _ X )f electrons around 40. Finally, basis set choice is shown to be important in the proposed procedure 
"■^^o include electron correlation effects in determining thermodynamic properties. With the optimum 
oice of basis set, the average percentage error for some 20 molecules is approximately 20% for heats 
^ pi formation. For molecular dissociation energies the average error is much smaller: ~ 0.4%. 

'^%eywords: Molecular equilibrium energies; Molecular dissociation energies; Molecular heats of 
f~]f nrmation. 
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1 Introduction 

In the last few years, much effort has been devoted to the development of theoretical techniques for 
the calculation of fundamental thermodynamical quantities related to molecular formation. These 
quantities include the molecular heat of formation, A.H , the variation of entropy, AS, and the 
variation of free energy, AG, as well as other physical quantities related to them, such as the molecular 
dissociation energy Do, and the equilibrium constants of the chemical reactions. 

The available theoretical techniques to extract such quantities are especially helpful whenever the 
experimental results, albeit accurate, are affected by uncertainties, due to the difficulty of obtaining 
pure samples of a given compound. Therefore, efScient and accurate theoretical methods are required 
to compute thermodynamical quantities. 

In particular, the Gx methods {x = 1,2,3) have been recently developed to this purpose U. 
These techniques are based on the combination of several ab initio molecular energy calculations, 
using various basis sets, and including correlation effects within the framework of M0ller-Plesset (MP) 
theory 0, at MP2 and MP4 levels of approximation. However, one of the main drawbacks of the 
Gx methods is their computational complexity. Indeed, even for the energy calculation of molecules 
of relatively small size, both computation time and memory occupancy are quite expensive As 
a consequence, the application of this class of methods is effectively limited by the system size. 
Therefore, it is of considerable interest to investigate alternative theoretical techniques to calculate 
useful thermodynamical quantities for the formation of molecules, which require less computational 
resources, while maintaining a good level of accuracy. 

In this context, Cremer et al. I5lbl have demonstrated that it is possible to extract the Schrodinger 
energy of a molecule, Es, i.e. the energy corresponding to the exact solution of the Schrodinger 
equation for a given molecule, from the experimental values of some observable quantities, such as 
the molecular heat of formation and the molecular vibrational frequencies. Inversely, any theoretical 
model able to closely reproduce the molecular Schrodinger energies, is also expected to yield molecular 
heats of formation AH close to the experimental ones. 

By definition, the difference between the molecular Schrodinger energy Es and the molecular 
energy obtained from a calculation at Hartree-Fock (HF) level is the 'experimental' correlation energy, 
E^^'^. For this reason, in Ref. [7] we developed a semi-empirical method to calculate a theoretical 
estimate iJ^hco ^^j^g bond-order correlation energy, or BOCE) to the correlation energy Ec^'^ based on 
the molecular bond-order matrix. Within the framework of the Lewis model ,8,, which describes the 
bond between atoms in a molecule as the sum of localized electron pairs, this 'theoretical' correlation 
energy £;*hco jg calculated as the sum of the correlation energies of each atom, and the contribution 
due to the binding energy of each bond. The latter is in turn evaluated in terms of the product 
between the molecular bond-order matrix and some empirical parameters (see also Ref. ^ for more 
details). Since the BOCE method relies on a relatively simple transformation of the matrix of the 
molecular orbital coefficients, it turns out that it is a rather fast and accurate technique, as compared 
to the Gx methods. 

The BOCE technique has been applied to calculate molecular properties, such as the molecular 
energy E^^'^°, the ionization potentials, and the electron affinities of closed and open shell molecules 
containing C, H, and O atoms ,7., and recently it has been extended to closed and open shell molecules 
containing Si, H, F, and CI .9|. In view of the previous satisfactory results for these molecular classes, 
in this paper we apply the BOCE method to the calculation of other molecular properties such as the 
molecular heat of formation and the molecular dissociation energies of several molecules containing 
C, H, and O. Since the bond-order matrix depends on the basis set used in the HF calculation, we 
compare the results within the G2 model with our BOCE results obtained using different basis sets. 



2 Method 

The starting point of the BOCE method is the experimental correlation energy, E"^'^. This is defined 
as the difference between the Schrodinger molecular energy Eg and the Hartrce-Fock energy, -EhFi 



E^^'P — Eg — i?HF- 



(1) 



In Eq. 0, -BjjF is the energy from an ideal Hartree-Fock calculation, requiring an infinite basis set. 
It should be emphasized that, in any practical calculation, -Bhf includes errors due to the particular 
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choice of basis set in the HF calculation. These errors can be made arbitrarily low, by choosing a 
sufficiently large basis set. 

In analogy to Eq. one can define a 'theoretical' correlation energy E^^°°, as the difference 
between the molecular energy iJ'hco obtained within a theoretical model at higher level than HF 
(below, the superscript 'theo' will alternatively refer to Gx or BOCE), and the HF molecular energy: 



^thoo ^ £,thoo _ ^^p^ (2) 
Following the molecular dissociation scheme, reported in Ref. 

AB — ► A + B + binding energy, (3) 

both experimental and theoretical correlation energies can be partitioned into sums of atomic and 
bond contributions, 

= E '^""'"{AB), (4a) 

A all AB 

E*J^<=° = J2eI^''°{A) + J2 E*J^''°{AB), (4b) 
A all AB 

where the sums run over all atoms and pairs of atoms in the molecule, respectively, and both exper- 
imental and theoretical atomic correlation energies are defined as in Eqs. and J^Ji respectively. 

As assumed in Refs. I7I9I . within the BOCE method (theo = BOCE) we identify the theoretical 
energy of a single atom with its Schrodinger energy, 

Es{A) = E^'^'='^{A). (5) 

Therefore, the atomic theoretical correlation energies will be equal to the experimental ones: 

E^°'^^{A) = E^^''P{A). (6) 

On the other hand, the bond contribution within the BOCE method is described as an analytical 
function of the bond order Pas between atoms A and B, as obtained within the HF calculation. 
Therefore, to lowest order in Pab, one has EE) 

£;BOCE(^^) ^^^^p^^_ (7) 

where a^B is a parameter which depends on the A — B bond. Summarizing, the theoretical molecular 
correlation energy within the BOCE method can be expressed as 



^BOCE 



■.J2e^^p{A)+ J2 ^abPab- (8) 

A all AB 



The bonding parameters oab in the second term of Eq. JSJ can be obtained from the experimental 
correlation energies of model molecules. In particular, having M different pairs of atoms in a molecule, 
each of them contained ra^ times, we can rewrite the second term of Eq. JSJ as 



y , iabPab = niaABPAB + n2aAcPAC + ■ ■ • + nMawzPwz, 
all AB 



(9) 
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where it has been assumed that the bond-order of each pair is approximately the same for each pair 
of the same kind in the molecule. In this sense, the parameter aAB is the binding correlation energy 
per bond A — B. 

The procedure starts with the calculation of the parameter uhh obtained from the experimental 
correlation energy of a H2 molecule, 

Er^(n2)-2ET^{K) .^^^ 
'^™ = 1^ ■ ^'"^ 

In a similar way, the parameter acH is obtained from the experimental correlation energy of the CH4 
molecule, and making use of the value for ohHi 

i?r''(CH4) - [4£;rP(H) + ET^iC) + eaHH^HH] 
acH = 



4PcH 

grP(CH4) - 6aHHPHH 
4PCH 



(11) 



Following an analogous procedure, the parameters acCi ''COi ^OOi ^-i^d aoH have been obtained 
starting from the experimental correlation energies of C2H6, (CH3)20, H2O2, and H2O, respectively 
(see Table 

In concluding this section, let us briefly recall the definition of the thermodynamical quantities 
which will be calculated in the next section, viz. the molecular dissociation energy and the molecular 
heat of formation. 



2.1 Definition of some relevant thermodynamic quantities 

The experimental molecular dissociation energy is defined as 

= J2^siA)-Es + EzPE, (12) 



where Eg (A) is the Schrodinger energy of atom A in the molecule. Eg is the total Schrodinger energy 
of the molecule, and Ezpe is the vibrational Zero Point Energy, which is calculated from the values 
of the molecular vibrational frequencies. In view of Eq. JSJ, within the BOCE approximation one 
has for the molecular dissociation energy: 

D^O^^=Y,Es{A)-E^°^^ + EzPE. (13) 

A 

Therefore, the difference between the experimental and calculated dissociation energies can be ex- 
pressed as 

^cxp _ ^BOCE ^ ^BOCE _ 

_ ^BOCE _ ^cxp 

= [E'^'^'^^iAB) - Ef"P{AB)]. (14) 

all AB 

In other words, only the difference between bonding BOCE and experimental energies contributes 
to D^'P - DBOCE 
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Finally, following Cremer et al. I5IKI . the experimental molecular heat of formation is related to 
Es by 




and therefore the BOCE molecular heat of formation is related to eBOCE 

^ (^^^^ ^.Hok)~{j: Es{A) - bBOCe + ^^^^ j 

= (H29SK-HoK)~Dr''''- (16) 

In the above equations, -ff298 K ~ ^0 K the variation of molecular enthalpy from 298 K to K, 
and can be evaluated according to standard thermodynamical formulas 1101 . Therefore, subtracting 
Eqs. I15i and I16i . one obtains 

^^cxp _ ^^BOCE ^ £,BOCE _ ^oxp ^^^^ 

We recall that the above derivation holds in the BOCE approximation, where Es = iJ^OCE^ £qj, 
atoms. In general, for other methods, such as the G2 method, one has Es + E'^'^, and Eqs. 031 and 
I14i become; 



D^^ = Y,E'^\A)-E^^ + EzPE, (18) 

A 

and 

DT" - or = Y.^Es{A) - E^\A)] + (Es - E^% (19) 



respectively, whereas Eq. 1171 holds also for the G2 method, i.e. 

^^exp _ ^^G2 ^ ^G2 _ ^cxp ^20) 



3 Results and discussion 

The procedures involved in all Gx calculations {x = 1, 2, 3) are quite similar, and have been described 
in the original works by Curtiss et al. 0. With respect to the calculation of the molecular heats of 
formation, the results of the Gx methods show that, in some classes of compounds, the G3 method 
is more accurate, while for other classes of compounds the G2 or the Gl methods seem to be more 
accurate than the G3. We have chosen to compare the molecular heats of formation calculated within 
the BOCE approximation with those obtained with the G2 method. 

In all calculations involving the BOCE method, we have used the basis sets 3-21G, 6-311G, 
6-31G**, 6-311G**, and 6-311GH — h**, and in the calculation of the bond-order matrix we have 
employed Lowdin's definition )11| for all the latter basis sets. Our results for the bond parameters aAB 
for the different chosen basis sets are reported in Tab.lTl while our results for the bonding correlation 
energies, E^'~"^^{AB) are reported in Tab. 1^ for all pairs AB between H, C, and O. It may be 
seen that the variations in the bond parameters from one basis set to another are closely correlated 
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with those in the bonding correlation energies for the same bond and between the same basis sets. 
This is in agreement with the procedure outlined in the previous section [see e.g. Eq. Illl l. in view 
of the minor dependence of aAB on the bond-order matrix element Pab than on the experimental 
molecular correlation energy _E™^. 



3.1 Molecular energies 

Tab. 13 lists the set of 23 molecules, containing C, O, and H, considered in this work. For these 
molecules, Tab. l4l reports the experimental values of the Schrodinger energies, the total molecular 
energy (including correlation) within the G2 model, E'-^^ , and the total molecular energy (again, 
including correlation) within the BOCE approximation, ijBOCE _ j^BOCEj^j^^^ The latter depends 
on the particular basis set chosen for the HF calculation, as described above. 

It may be seen that the difference \E'^'^ — Eg\ ranges between fs! 140 and 535 mhartree, whereas 
j^BOCE _ jg a,lways below fs! 25 mhartree, for all the basis sets considered. Our first conclusion 
is therefore that, with respect to the Schrodinger result, the BOCE approximation is more accurate 
than the G2 method for the 23 molecules in our data set. While it is tempting to assume that 
this conclusion will remain true for the entire G2 data set, this is not decisive based solely on our 
present study. In particular, from Tab. |2 the best BOCE results correspond to the 6-311G basis 
set. This is due to a competition of various effects. First of all, the HF contribution to the total 
energy ijBOCE _ ^BOCE ^ _Ejjp clearly decreases on increasing the size of the basis set, owing 
to the variational nature of the HF calculation. Therefore, one invariably fi nds t hat -Ehf decreases 
along the series 3-21G — ► 6-311GH — h**. On the other hand, following Eq. 14bl . and decomposing 
the correlation energy £;BC)CE jj^^-^ atomic and binding contributions, one finds that the atomic 
(respectively, binding) contributions for the various basis sets are minimum (respectively, maximum) 
for the 6-311G basis set. This is probably due to the nature of the 6-311G basis set (and of the 6-3xx 
basis set family, in general), usually employed to optimize the energies of the single atoms. 

As a result of the compensation of these various contributions, the use of the 6-311G basis set 
within the BOCE method yields the best agreement with the Schrodinger energy. This is especially 
important, in view of the dependence of the dissociation energies Do on Es. 

In concluding this subsection on the molecular energies, let us discuss the dependence of the bond 
order for a given bond on the basis set employed. The average values of the bond order Pab for 
the bonds C— H, C— C, O-H, C-O appearing in the molecules listed in Tab. ISl have been reported 
in Tab. 13 for the various basis sets. (The bond O-O has not been included, since there is only one 
molecule in Tab. l3l involving such bond.) One can see that, in bonds involving H and C atoms, the 
bond order variation along the series of basis sets considered here is below 4%, whereas in bonds 
involving the oxygen atom, the bond order variat ion ranges from ^ 18% (C— O) to Ri 27% (O— H). 
This is in agreement with the study of Sannigrahi [121 . where a more pronounced dependence of the 
bond-order and of the valency on the basis set is found for all compounds containing electronegative 
centers (such as O, F, CI). For calculations concerning these compounds, Sannigrahi 1121 therefore 
recommends the use of double-zeta basis sets, and of polarization and diffuse functions in the case 
of highly ionic and/or negatively charged species. 



3.2 Dissociation energy and molecular heat of formation 

In Tab.|S]we report the experimental and theoretical molecular dissociation energies Dq for the set 
of molecules listed in Tab. l3l In particular. Tab. l6l compares the theoretical values obtained within 
the G2 method (second column) and the BOCE approximation, for the five basis sets considered 
above. Analogously, Tab. Q reports the experimental and theoretical molecular heats of formation 
AH for the same molecules, models, and basis sets. 

In the case of the theoretical dissociation energy calculated within the G2 method, , we find 
an average percentage difference with respect to the experimental value of fa 0.2 % (corresponding to 
fa 1.5 kcal). For the same quantity calculated within the BOCE approximation, we find an average 
percentage difference with respect to the experimental value of Ri 0.4 — 0.8 % (corresponding to 
Si 2.4 — 5.9 kcal), depending on the basis set employed (Tab. 1^. Analogously, for the molecular 
heats of formation (Tab. 0, on the average we find |AHC=2 _ A//<=''P| S3 8.2 % (corresponding to 
K3 1.5 kcal), and lA/fBO^E _ A//™p| 21.4 - 58.5 % (corresponding to « 2.4 - 5.9 kcal). (Such 
larger values of the average percentage errors are justified by the smaller values of the experimental 
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molecular heats of formation, than the experimental molecular dissociation energies.) In both cases, 
the best agreement within the BOCE approximation is obtained for the 6-311G basis set, as it was 
the case for the molecular energies. 

Therefore, we may conclude that, while the G2 method yields in general more accurate estimates 
of the molecular dissociation energies and the heats of formation, the BOCE approximation affords 
theoretical estimates of the above quantities of comparable accuracy, but now requiring much less 
computational effort. Within the BOCE approximation, moreover, the best agreement with the ex- 
perimental results is obtained for the 6-311G basis set. As observed for the molecular energies, this 
is a result of the competition of the various atomic and binding contributions to the correlation ener- 
gies (and therefore to the other thermodynamical quantities considered in this work). In particular, 
within the BOCE approximation, the binding contribution to the total correlation energy is relatively 
larger than the atomic contribution, and it can be therefore directly related to the bond-order value. 

While it is agreed that the use of modest basis sets reduce the computing time required to obtain 
chemical accuracy relative to the Gx methods, it is appropriate, at this point, to expand on our 
reasoning as to why the 6-311G basis set gives better results than the somewhat larger basis sets 
we have also worked with. For such, still quite restricted, basis sets, there is rather non-uniform 
'convergence'. We offer, as an explanation of such non-uniformity, that the better results obtained 
with the 6-311G basis set rest on a somewhat delicate interplay between diffuse s-like (-I-) and 
polarization d-like (*) functions on the one hand, and functions corresponding to valence electrons 
on the other. As a consequence, in the larger basis sets considered here, the total bond-order of each 
atom is overestimated with respect to the 6-311G basis set. The latter point about bond-order can be 
illustrated with reference to Table l5l If we neglect the 3-21G basis set recorded there, then with just 
two exceptions (the C-C bond with the 6-311G** basis set and the C— H bond with the 6-311G-I — h** 
basis set), the bond orders are smallest for the 6-311G set. Specifically then, the advantage among 
these small basis sets lies with the 6-31 IG basis set in that its diffuse and polarization functions 
increase, in the final bond-order matrix, the total number of electrons in each atom contributing to 
the bond. 

A finer analysis of the various terms contributing to the molecular dissociation energies and heats 
of formation reveals that the very good agreement of the G2 method with the experimental results is 
due to a partial compensation of the two competing contributions for D'^^ — D'^^ and AH™P — Ai?'^^ 
in the right hand sides of Eqs. 1191 and 1201 . respectively. Indeed, one finds E"^^ - Es ~ 74.356 - 
229.413 kcal, while EaI-^'^^C^) " ^siM ~ 74.169 - 222.507 kcal, for the various molecules of 
Tab. 13 with differences between the two quantities ranging between —3.0 kcal and 6.9 kcal (Fig. Q. 
Therefore, in some way, within the G2 method the errors on the molecular part are compensated by 
the atomic part, i. e. the binding energy is very close to the Schrodinger value. 

On the contrary, within the BOCE method [cf. Eqs. 1141 and II17I 1 . the difference between exper- 
imental and calculated molecular dissociation energies and heats of formation depends only on the 
difference between the Schrodinger and the BOCE molecular energy, given the assumption that the 
BOCE atomic energies are equal to the Schrodinger counterparts, Eq. JSJ. Therefore, even though 
the BOCE molecular energy is very close to the Schrodinger value, such a small difference (in kcal) 
cannot be compensated by the atomic contribution. 



4 Conclusions and directions for future work 

We have compared and contrasted the Gx method (especially with x = 2) and the BOCE approxi- 
mation for estimating (a) molecular energies, (b) molecular dissociation energies, and (c) molecular 
heats of formation for some 20 molecules containing C, H, and O atoms. Although the G2 method 
usually yields better agreement with the experimental results, we find that the accuracy of the the- 
oretical estimates of the above observable quantities within the BOCE approximation is comparable 
to that of the G2 calculations, but now with a remarkable saving in terms of computational com- 
plexity. With respect to previous studies I7I9I . all BOCE parameters have been calculated on single 
model molecules, rather than averaging over several molecules. Therefore, each BOCE parameter 
OAB refiects the nature of the particular bond A — B, irrespective of its chemical surroundings. In 
other words, we find that the BOCE method is rather robust with respect to the inclusion of further 
than binary correlation terms. Moreover, we have extensively analyzed the dependence of the BOCE 
results on the basis set chosen for the underlying HE calculation. 

Owing to the much reduced computational effort required in a BOCE calculation, we plan to 
apply the latter method to the evaluation of molecular properties of larger molecules, where Gx 
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methods are expected to require presently prohibitive resources. It is also intended, in the light of 
the excellent results for molecular dissociation energies, to study for such a set of larger molecules 
the relation between the findings of the BOCE approximation and Teller's theorem [T21, which states 
that molecules do not bind in a fully local density approximation (LDA), i.e. including kinetic energy 
density in LDA. Mucci and March 14, proposed to make a merit out of Teller's theorem, by relating 
dissociation energy to molecular electron density gradients llSllbl . However, electron correlation in 
the separated atoms seems important also in respect to heats of formation [see Eqs. I19i and 1201 
above], and this will also require further exploration. 
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Tabic 1. Bond-order parameters a^B in Eq. (7J (in mhartree) for all possible bonds A — B between C, O, 
and H. Their values depend on the HF reference energy, which depends in turn on the particular choice of 
basis set. 



A- 


- B 


Model molecule 


3-21G 


6-311G 


6-31G** 


6-311G** 


6-311G++** 


H 


H 


H2 


43.86 


46.00 


39.55 


41.55 


41.55 


C- 


-H 


CH4 


39.79 


41.51 


35.12 


34.85 


35.59 


C- 


-C 


C2H6 


39.77 


39.60 


31.70 


30.68 


30.06 





-H 


H2O 


88.43 


84.63 


60.54 


51.20 


51.71 


C- 


-0 


(CH3)20 


67.85 


68.14 


49.11 


43.91 


41.11 


0- 


-0 


H202 


83.33 


89.52 


72.04 


69.92 


66.62 



Table 2. Bond correlation energies, Ec{AB), in hartree, for all possible bonds A — B between C, O, and H. Their values depend on the HF reference energy, 
which depends in turn on the particular choice of basis set. 



A- 


-B 


Model molecule 


3-21G 


6-311G 


6-31G** 


6-311G** 


6-311G++** 


H- 


-H 


H2 


-1.17442 


-1.17442 


-1.17442 


-1.17442 


-0.17442 


C- 


-H 


CH4 


-40.51353 


-40.51353 


-40.51353 


-40.51353 


-39.51353 


C- 


-C 


C2H6 


-79.82261 


-79.82261 


-79.82261 


-79.82261 


-78.82261 


0- 


-H 


H2O 


-76.43076 


-76.43076 


-76.43076 


-76.43076 


-75.43076 


c- 


-0 


(CH3)20 


-151.54948 


-151.54948 


-151.54948 


-151.54948 


-150.54948 


0- 


-0 


H202 


-155.02013 


-155.02013 


-155.02013 


-155.02013 


-154.02013 
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Table 3. List of the 23 molecules, containing C, O, H, considered in this work. 



1 

1 


benzene 




9 


eijuyieiie 




3 




C2H2 


A 


iUi ilidlU-Clij U-C 


TTPHO 


D 


VKl /!!\"r r^T T" 1 o 1 /~\ h» /~\ 1 

ineijiiyi dicunui 




Q 


JVC UCllC 




7 


i^Q-prAi^tn (HirwiHp 






o 


(dccid-iciciiycie 


^n3V_>jrLVj 


Q 


CLllCIlOl 


^n2^xlv>'xl 


1 n 


formic 3,ci(i 


I 1 V V / w 1 1 


1 1 

± X 






1 9 


cycloprop3-ne 


(,*-'ll2j3 




1 , Zi pi (jpdLiieiie 


vyn2 v_^v^JrL2 


14 


fur an 


C4H4O 


15 


cyclohexane 


(CH2)6 


16 


glyoxal 


HCOCOH 


17 


1,3 butadiene 


CH2CHCHCH2 


18 


acetone 


CH3COCH3 


19 


acetic acid 


CH3COOH 


20 


propene 


CH3CHCH2 


21 


2-butene (E) 


t-CH3CHCHCH3 


22 


2-butene (Z) 


C-CH3CHCHCH3 


23 


2-methyl 1-propene 


C(CH3)2CH2 



Tabic 4. Total molecular energies (in hartroo), £;B0CE ^ ^bocb _|_ ^^^^ including correlation, for the 23 molecules listed in Tab. The various columns 
refer to the different basis sets being considered in the HF calculation. The second and third columns are the total molecular Schrodinger energy, Es, and 
, respectively. 





Es 
















3-2 IG 


6-311G 


6-31G** 


6-311G** 


6-311G++" 


1 
1 


—232 


24219 


—231 


78054 


—232 


23115 


—232 


24129 


—232 


22161 


—232 


22495 


—232 


21644 


9 


—78 


58567 


—78 


41594 


—78 


58344 


—78 


58651 


—78 


58084 


—78 


58302 


—78 


58157 


o 
o 


—77 


33430 


—77 


18574 


—77 


33024 


—77 


33251 


—77 


31958 


—77 


32522 


—77 


32114 


A 


— 114 


50344 


— 114 


O O O A O 

33892 


— 114 


49660 


— 114 


50112 


— 114 


50278 


— 114 


50176 


— 114 


A r\ A AA 

49490 





— 115 


72137 


— 115 


53489 


— 115 


72411 


— 115 


72416 


— 115 


72340 


— 115 


72350 


— 115 


TO A AO 

72202 


D 


— 152 


6 


— 152 


36911 


— 152 


60331 


— 152 


59168 


— 152 


59566 


— 152 


59897 


— 152 


58538 


7 


— 188 


58504 


— 188 


36132 


— 188 


58684 


-188 


57255 


-188 


60229 


— 188 


60371 


-188 


58161 


Q 
O 


— 153 


O O A O /t 

82924 


— 153 


57685 


— 153 


82178 


— 153 


82718 


— 153 


82682 


— 153 


82746 


— 153 


O O A ^ A 

82040 


Q 


— 153 


81314 


— 153 


55770 


— 153 


81266 


— 153 


81354 


— 153 


O A O O O 

80838 


— 153 


81022 


— 153 


80530 


1 {^ 


— 189 


76370 


— 189 


51649 


— 189 


77233 


-189 


76925 


-189 


77838 


— 189 


77720 


-189 


76468 


1 1 
1 ± 


— 113 


31751 


— 113 


17750 


— 113 


31573 


— 113 


31644 


— 113 


O A A A O 

32008 


— 113 


32148 


— 113 


31558 


1 9 


— 117 


89065 


— 117 


63115 


— 117 


87292 


— 117 


O O A A O 

88093 


— 117 


O O A O A 

88280 


— 117 


88375 


— 117 


88177 


1 


— 116 


65212 


— 116 


41784 


1 1 fi 

— 1 ±D 




lift 




lift 


uoyuo 


— ± ID 




— 1 ID 


ooooy 


1 A 


— zou 


uiyou 


99Q 


UOZUl 


-230 


01332 


-230 


01476 


-230 


00899 


-230 


00925 


-229 


99859 


15 


-235 


87496 


-235 


33937 


-235 


85806 


-235 


86391 


-235 


86104 


-235 


86280 


-235 


86424 


16 


-227 


81908 


-227 


51024 


-227 


81254 


-227 


82108 


-227 


82216 


-227 


82183 


-227 


80646 


17 


-155 


98889 


-155 


66427 


-155 


97929 


-155 


98644 


-155 


97466 


-155 


97879 


-155 


97500 


18 


-193 


15209 


-192 


81368 


-193 


14508 


-193 


15078 


-193 


14894 


-193 


14943 


-193 


14242 


19 


-229 


08873 


-228 


75391 


-229 


09387 


-229 


09385 


-229 


09930 


-229 


09921 


-229 


08686 


20 


-117 


90274 


-117 


64509 


-117 


89730 


-117 


90190 


-117 


89541 


-117 


89816 


-117 


89658 


21 


-157 


21942 


-156 


87484 


-157 


20633 


-157 


21660 


-157 


21007 


-157 


21245 


-157 


18493 


22 


-157 


21848 


-156 


87266 


-157 


20984 


-157 


21412 


-157 


20760 


-157 


21001 


-157 


20536 


23 


-157 


22188 


-156 


87670 


-157 


21280 


-157 


21716 


-157 


21044 


-157 


21311 


-157 


18634 
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Table 5. Dependence of the bond-order Pab for a given bond A — B on the basis set employed. Reported 
are the average bond-orders Pab for the bonds more frequently appearing in the molecules listed in Tab.lsl 





3-21G 


6-311G 


6-31G** 


6-311G** 


6-311G++** 


c- 


-H 


0.94649 


0.93456 


0.94714 


0.95081 


0.91798 


c- 


-C 


1.05140 


1.06971 


1.07062 


1.05695 


1.09041 


0- 


-H 


0.92589 


0.91894 


1.00533 


1.16700 


1.12736 


c- 


-0 


1.32735 


1.15963 


1.24610 


1.32735 


1.37367 



Tabic 6. Experimental and calculated dissociation energies Dq (in kcal) for the 23 molecules listed in 
Tab. 1^ Boldface values refer to the method yielding the best agreement with the experimental result, 
for each molecule. The two last lines refer to the average absolute difference and the average percentage 
difference, respectively, between the experimental and calculated Dq. 







nG2 
Dq 










3-21G 


6-311G 


6-31G" 


6-311G" 


6-311G++" 


1 


1305.5 


1301.8 


1298.6 


1305.0 


1292.6 


1294.7 


1289.4 


2 


531.9 


531.7 


530.5 


532.4 


528.8 


530.2 


529.3 


3 


388.8 


387.2 


386.2 


387.6 


379.5 


383.1 


380.5 


4 


359.0 


359.3 


354.7 


357.6 


358.6 


358.0 


353.7 


5 


480.9 


482.3 


482.6 


482.6 


482.1 


482.2 


481.3 


6 


512.8 


513.6 


516.6 


509.3 


511.8 


513.8 


505.3 


7 


381.9 


384.6 


383.1 


374.1 


392.8 


393.7 


379.8 


8 


643.7 


643.9 


639.0 


642.4 


642.2 


642.6 


638.2 


9 


633.6 


631.9 


633.3 


633.9 


630.6 


631.8 


628.7 


10 


480.0 


482.0 


485.4 


483.4 


489.2 


488.4 


480.6 


11 


256.2 


258.0 


255.1 


255.5 


257.8 


258.7 


255.0 


12 


802.9 


802.1 


791.7 


796.8 


797.9 


798.5 


797.3 


13 


669.1 


668.2 


665.1 


667.4 


661.2 


663.7 


660.7 


14 


950.5 


949.6 


946.4 


947.3 


943.7 


943.9 


937.2 


15 


1659.4 


1652.5 


1648.8 


1652.5 


1650.7 


1651.8 


1652.7 


16 


610.4 


613.4 


606.3 


611.7 


612.4 


612.2 


602.5 


17 


960.1 


958.2 


954.1 


958.6 


951.2 


953.8 


951.4 


18 


927.0 


927.9 


922.6 


926.2 


925.0 


925.3 


920.9 


19 


764.7 


766.3 


767.9 


767.9 


771.3 


771.3 


763.5 


20 


811.2 


810.8 


807.8 


810.6 


806.6 


808.3 


807.3 


21 


1090.5 


1090.3 


1082.3 


1088.7 


1084.9 


1086.2 


1084.9 


22 


1089.8 


1089.0 


1084.4 


1087.0 


1083.2 


1084.7 


1083.2 


23 


1092.2 


1091.5 


1086.5 


1089.2 


1085.1 


1086.7 


1085.9 


Avg 


|Ai?ol 


1.5 


4.5 


2.4 


5.5 


4.6 


5.9 


Avg 


lAi^ol (%) 


0.2 


0.6 


0.4 


0.8 


0.7 


0.8 



Table 7. Experimental and calculated molecular heats of formation at 298 K, Ai? (in kcal), for the 23 molecules listed in Tab.[^ Boldface values refer to 
the method yielding the best agreement with the experimental result, for each molecule. The two last lines refer to the average absolute difference and the 
average percentage difference, respectively, between the experimental and calculated /S.H . 
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6-311G 


6-31G** 


6-311G** 


6-311G++** 
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1 Q 8 




23.6 


26 8 


20.4 


32 7 

cJZ< . J 


30 6 

KjyJ .\j 


36 


2 


12.5 




12.8 


13.9 


12.0 


15.6 


14.2 


15.1 


3 


54.2 




55.8 


56.7 


55.3 


63.4 


59.9 


62.5 


4 


— 27 7 




-28.0 


—23.4 


—26.2 


—27 3 

Zj i .O 


—26 7 

ZU . 1 


—22.3 


5 


—48 n 

4:0. U 




-49.4 


—49 8 


—49 8 


—49 3 


—49.4 


—48.4 


6 


— lis 




-12.2 


— 15.2 


— 7 9 


— 10.4 


— 12.5 


—4 


7 


— 94 1 




-95.4 


—95.2 


— 86 2 


— 1 04 9 




— 91 9 


8 


—40 8 




41.0 


— 36 1 


— 39 S 


— 39 3 


— 39 7 


— 3^ 3 


9 


—30.6 




-29.0 


—30.3 


—30.8 


—27.6 


—28.8 


—25.7 


10 


— 90 5 




-92.6 


—95 9 


—94 n 


—99 7 


—99 


— Q1 1 


11 


—26.4 




-28.2 


— 25 3 




— 28 


— 28 9 

ZO . J 


—25.2 


12 


1 2 7 




13.5 


23 9 


1 8 8 


1 7 7 


1 7 1 


1 8 3 

iO .O 


13 


4S 5 




46.4 


49.4 


47.2 


53 4 


•lO 9 


^3 8 


14 


— 8 3 




-7.4 


—4.2 


— 5 1 


— 1 5 


— 1 7 


S 


15 


-29.4 




-22.5 


-18.8 


-22.5 


-20.7 


-21.8 


-22.7 


16 


-50.7 




-53.6 


-46.6 


51.9 


-52.6 


-52.4 


-42.7 


17 


26.0 




28.0 


32.0 


27.5 


34.9 


32.3 


34.7 


18 


-52.2 




53.0 


-47.8 


-51.4 


-50.2 


-50.6 


-46.2 


19 


-103.44 




-105.1 


-106.7 


-106.7 


-110.1 


-110.0 


-102.3 


20 


4.9 




5.3 


8.3 


5.4 


9.5 


7.8 


8.7 


21 


-2.6 




2.4 


1.8 


-0.8 


3.0 


1.7 


3.0 


22 


-1.8 




1.0 


3.6 


0.9 


4.8 


3.4 


4.7 


23 


-4.3 




5.0 


1.4 


-1.2 


2.9 


1.2 


2.1 


Avg 




- AiJ^^'P 


1.5 


4.3 


2.4 


5.5 


4.6 


5.9 


Avg 




- AiJ^'^P] (%) 


8.2 


44.3 


21.4 


52.9 


41.8 


58.5 
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Figure 1. Atomic vs molecular contributions to the difference between experimental and 
theoretical dissociation energy, Eq. 1191 . within the G2 method, for the 23 molecules listed i 
Tab. El The dashed line is a guide to the eye. 



